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Abstract 



We study the thermodynamic hmit of the particle-hole form factors of the 
XXZ Heisenberg chain in the massless regime. We show that, in this limit, 
such form factors decrease as an explicitly computed power-law in the system- 
size. Moreover, the corresponding amplitudes can be obtained as a product 
of a "smooth" and a "discrete" part: the former depends continuously on the 
rapidities of the particles and holes, whereas the latter has an additional explicit 
dependence on the set of integer numbers that label each excited state in the 
associated logarithmic Bethe equations. We also show that special form factors 
corresponding to zero-energy excitations lying on the Fermi surface decrease as 
a power-law in the system size with the same critical exponents as in the long- 
distance asymptotic behavior of the related two-point correlation functions. 
The methods we develop in this article are rather general and can be applied 
to other massless integrable models associated to the six-vertex i?-matrix and 
having determinant representations for their form factors. 



1 Introduction 

This article is a continuation of our study of form factors in massless quantum integrable 
models. We have recently flj investigated the thermodynamic limit of a special, Umklapp- 
type, form factor of the XXZ Heisenberg spin chain OEIIHEIEIITIIHIIH] starting from its 
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determinant representation [lUj in the algebraic Bethe Ansatz framework [IH \H[ [T^ . We now 
apply this method to general, particle-hole type, form factors. The main goal of our analysis is 
to pave a way for the calculation of the asymptotic behavior of correlation functions in massless 
integrable models through their form factor expansion. Our results should also be useful in 
the numerical computation of dynamical structure factors through their (exact) form factor 
expressions starting from finite size systems [14^ [T5| [TBI \T7\ [T8] as they allow a precise a priori 
control of the dependence of each form factor in terms of the system size. 

For definiteness, we chose to focus on a particular model, the XXZ spin- 1/2 Heisenberg 
chain [2] in the massless regime and in an external magnetic field h > for which determinant 
representations of general form factors were obtained in [lOj. However, up to minor modifica- 
tions specific to the choice of the model, our method and results apply as well to other massless, 
algebraic Bethe Ansatz-solvable models with known determinant representations of their form 
factors like the non-linear Schrodinger model [19^ [20] or the higher spin XXX [21] and XXZ 
p2j chains. 

The Hamiltonian of the XXZ chain is given by, 

M , M 

H = J2 i^k^k+i + ^M+i + - 1)) - 2 E (1-1) 

fc=l k=l 

where cr^'^'^ are the spin operators (Pauli matrices) acting on the /c*^ site of the chain. The 
main purpose of this article is to study the asymptotic behavior of the form factors in the 
thermodynamic limit M — )• oo of the chain and in the massless regime |A[ < 1. We will take 
this limit starting from a chain of finite even size and subject to periodic boundary conditions. 

Considering the finite chain allows us to define rigorously the form factors of local spin 
operators as the normalized matrix elements 



where | V' ) and | tp') are two eigenstates of the Hamiltonian (jl.ip . Such matrix elements can be 
computed in a systematic way by using the solution to the quantum inverse scattering problem 
[TUl [23] together with the explicit formulae for scalar products [211 [ISl [TO] in the algebraic 
Bethe Ansatz framework pTl [T2l [T3l [26] . 

Using the closure relation, one can express any zero temperature two-point correlation 
function of spin operators as a sum over the form factors of the corresponding local operators, 

Here | ij^g ) denotes the ground state of the Hamiltonian (|l.ip . and the summation is taken with 
respect to all the eigenstates | -0') of (II. 1|) . 

One possible way of dealing with ()1.3p was proposed in [271 128j . where the sum over the 
form factors for the finite chain was recast into a multiple contour integral, the so-called master 
equation. This method then leads to various multiple integral representations for the correlation 
functions in the thermodynamic limit |27[ [28l [29] that are in fact non-trivial summations of 
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previously obtained multiple integral representations for the correlation functions [60\ [511 [5^ 
[33 | 1341 [35] . Moreover, in [36| we have shown how to derive long-distance asymptotic behavior 
of certain two-point correlation functions using such master equation representation. It agrees 
with predictions arising from the Luttinger-liquid and Conformal Field Theories approaches 
[371 EEl EHl Sol SH SSI Hal E]. it is the very remarkable structure of the results obtained in 
[36j (the amplitudes of the power law decrease of the long distance behavior of the two-point 
correlation function are related to the modulus squared of some properly normalized form 
factors) that strongly suggested that another possible way to analyze the asymptotic behavior 
of the two-point correlation functions would be to take the thermodynamic limit directly in 
p.3|) . It seems indeed quite natural that, in this limit, the main contribution to the sum (|1.3|) 
should stem from excitations above the ground state having a finite energy. Such excited states 
can be characterized in terms of particles, holes and/or string configurations. In this framework, 
the sum over the complete set of states in ()1.3p becomes a sum over all the possible particle- 
hole-string type excitations above the ground state which, after some possible regularizations, 
should be replaceable by integrals in the thermodynamic limit. Another important motivation 
to learn how to deal directly with the form factor expansion of the correlation functions concerns 
the time dependent case, as in that situation the asymptotic analysis through the dynamical 
master equation [28] poses yet unsolved problems. The first step to carry out this program 
is therefore to determine the leading asymptotic behavior of the form factors of the model in 
the thermodynamic limit. It is the purpose of the present article to solve this problem for 
the so-called particle/hole type form factors (corresponding to diffusion states). The case of 
states associated with complex solutions of the Bethe equations such as string type solutions 
(corresponding to bound states) requires the use of additional techniques and will be considered 
in a separate publication. 

In [I], we have shown that, in the M oo limit, a special Umklapp-type form factor of the 
cj^ operator decreases as some negative power of the system size M. The striking feature of this 
analysis is that this power coincides with one of the critical exponents appearing in the long- 
distance asymptotic behavior of the zero temperature two-point correlation function (<7fcr^^^). 
In the present article, we will show that general form factors exhibit similar properties: they 
all decay as some negative power of the size M; moreover, for special types of form factors 
corresponding, in the thermodynamic limit, to particle-hole excitations lying on the Fermi 
surface, these exponents coincide exactly with the critical exponents describing the power- 
law decay of the long-distance asymptotic behavior of the corresponding two-point correlation 
functions. 

More precisely, in the thermodynamic limit, the product of two form factors appearing in 
()1.3p can be presented in the following form 

• = e'''-''^"'^ Sss' Vss'. (1.4) 

All the dependence on the lattice distance m — m' is contained in an obvious phase factoi0. 
There Vex is the excitation momentum. The finite part of the form factors product can be 
separated into two parts Sgs' and Vgs'- We call them smooth and discrete parts respectively. 

^Some form factors (see e.g. ((311}) may also have an additional factor (— l)™^™ , which can be removed by 
the re-definition of the Hamiltonian 
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We show that the smooth part Sss' depends continuously on the rapidities of the particles and 
holes. On the contrary, it is only when the particles (holes) rapidities are separated from the 
Fermi boundary that the discrete part Vgs' does also depend smoothly on these quantities. As 
soon as the latter approach the Fermi surface, the discrete structure of the form factors reveals 
itself in Vgs'- This means that a microscopic (of order 1/M) deviation of a particle (hole)'s 
rapidity leads to a macroscopic change of Vss'- We believe that one of the consequences of such 
a structure is the particular role played by the Fermi boundary in bosonization techniques. 

The article is organized as follows. In section [2] we give a description of the low- lying excited 
states of the model. Following the ideas of [1], we consider excitations above the ground state 
of the twisted transfer matrix. This enables us to formulate precisely our results in section O 
i.e. the structure of the large size asymptotic behavior of general particle-hole form factors in 
the XXZ spin-1/2 chain. The remaining part of the article is devoted to the proof of this result. 
In section m we show how form factors can be reduced to scalar products of a special type that 
we can represent in terms of determinants for chains of finite length M. We then calculate 
the thermodynamic limit of the form factors of the operator in section [5l In particular, we 
show that if the rapidities of the particles and holes collapse on the Fermi surface then the 
form factors have a discrete structure. In section [6] we obtain analogous results for the form 
factors of the operator^. We would like to stress that several methods and results used 
in the present article are borrowed from our previous work |Tj. We conclude by discussing the 
possibility of applying our results to the calculation of the long-distance asymptotic behavior of 
the two-point functions at zero temperature. Several auxiliary Lemmas and proofs of technical 
character are presented in the two appendices. 

2 Space of states of the model 

Since our ultimate goal is to analyze the two-point correlation functions at zero temperature 
through (|1.3p . the form factors we consider here correspond to matrix elements of a local 
operator between the ground state and an excited eigenstate of the Hamiltonian (II. ip . In this 
section we recall how such states can be described in the large M limit. 

In the algebraic Bethe Ansatz framework [11^ [121 [T3] , the space of states of the model (in 
finite volume) is constructed by means of the Yang-Baxter algebra realized by the operator 
entries of the monodromy matrix. In the XXZ case, and more generally for the class of models 
associated with the six-vertex R-matrix, such monodromy matrices take the form 



where A, B, C, D are quantum operators depending on some spectral parameter A. It satisfies 
Yang-Baxter commutation relations given by the six-vertex R-matrix. In this framework, the 
eigenstates of the Hamiltonian (jl.ip coincide with the ones of the transfer matrix T(A) = 
tr T(A) = ^(A) + D(A), and can be parameterized as = |'!/'(Ai, . . . , Aat)), = 0, 1, . . . , M/2, 

^Recall that a± = i(a^ ± io"^) 
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in terms of a set {Aj}^^ of solutions to the logarithmic Bethe equations 



N 



MpoiXj) - ^{Xj - Xk) = 2Trnj, j = l,...,N. (2.2) 



k=l 



Here the functions po{X) and ^(X) are the bare momentum and phase, 

where < ^ < vr and cos(" = A. The numbers nj, —M/2 < rij < M/2, are integers (for 
odd) or half- integers (for even). 

2.1 Thermodynamic limit of the ground state 

The Bethe roots Xj, j = 1, . . . ,N, describing the ground state correspond to the solution of the 
system of logarithmic Bethe Ansatz equations (j2.2p with the special choice rij = j — {N + l)/2, 
where A^ is the number of down spins in the ground statfQ [i |13 El E]. Thus, for given M 
and h, these parameters Xj are fixed quantities. In order to describe how they behave in the 
thermodynamic limit {N,M — > oo such that N/M tends to some fixed average density D), we 
introduce the ground state counting function, 

N 

«.) = i^p„M-^j:«(.-A,) + i^, (2.4) 

k=l 

which is build so that C(Aj) = j /M for j = 1, . . . , A^. This function defines the discrete density 
of the minimal energy state in the A^ sector by p{ijj) = C'i^)- It is possible to argue that, in 
the thermodynamic limit, the parameters Xj fill densely a finite interval [— of the real axis 
(the Fermi zone), and that the discrete density goes to a smooth function p{uj) that solves the 
following integral equation 

<? 

P(A) + ^jK{X-p) Pip) dp = ^Po(A), with K{X) = t9'(A). (2.5) 

The value of the endpoint q of the Fermi zone is fixed by requiring that f^^piX) = D. Note 
that, since limM-i-oc? (Ai) = 0, the thermodynamic limit .^(A) of the counting function is the 
antiderivative of p{X) that vanishes at —q. Therefore ^(A) can be expressed in terms of the 
dressed momentum p(A), the Fermi momentum being defined hy kp = p{q) 

A 

p(A) = 27tJ p{p) dp, e(A) = b(A) + piq)]/27T, ttD = p{q) = kp. (2.6) 




^The number TV depends on the magnetization of the ground state, which is fixed by the overaU magnetic 
field h. Hereafter, we call A'' sector the subspace of the space of states having spins down. 
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Several other properties of the ground state can be described through the solutions of other 
linear integral equations. We introduce two such functions that will play an important role in 
our analysis, namely the dressed charge Z{X) and the dressed phase (p{X,i'), which satisfy 

Z{X) + ^ y K(A - /.) d/x = 1, (2.7) 

and 

q 

u) + ^jK{X- /x) (/>(/i, u) dfi = ^'&{X - v). (2.8) 
"<? 

In fact, these two functions are not independent. Using that K{X) = '&'{X), one easily obtains 
that 

Z{X) = l + <P{X,q)-<P{X,-q). (2.9) 
Another non-trivial relationship between Z and (J) has been established in [46^ I47j: 

l + ^{q,q)-cl){-q,q) = Z'\q). (2.10) 



2.2 Thermodynamic limit of the excited states 

To describe the other eigenstates of the transfer matrix T(A), it is convenient for further 
purposes to consider instead the twisted transfer matrix 7^(A) = A{X) + kD{X), where k 
is some complex (twist) parameter. Hereby the excited states of the Hamiltonian (jl.ip are 
obtained as the k — )• 1 limit of the eigenstates of the twisted transfer matrix Tk (A) [27^ [T] . 

In complete analogy with the standard algebraic Bethe Ansatz considerations, the eigen- 
states IV^kIIa*})) of Xt(A) can be parameterized by sets of solutions of the twisted Bethe equa- 
tions 

/ N -\- 1\ 

Mpo(/i^,)-5Z^K-^4) = 2^(^i-^^)+2vra, j = l,...,N,. (2.11) 

fc=i ^ ^ 

Such states |V'k({/^})) depend on A^^^; parameters {/^^j.}^]^, where is the number of spins 
down and we have set k = e^'^*". In this article we consider a to be a real number. This 
restriction $5(0) = is not crucial, but convenient, since then (^^({/w})! = ("1)^" IV'k({/^}))^- 
Remark 2.1. One could of course from the very beginning set a = and deal directly with the 
standard excited states of the XXZ Hamiltonian (II. ip . However, we chose to introduce this 
extra parameter and keep it throughout our computations since, as we will show later on, an 
important class of form factors can be obtained by a mere shift of a by some integer value. 

In the following, we will consider two cases of interest: A^ = A for form factors of , and 
Ak = a + 1 for the form factor^. It follows from [8] that the set of real roots of the twisted 
Bethe equations is completely defined by the choice of the set of integers ij in the r.h.s. of 
()2.1ip . Therefore, we have labeled the roots of (12. lip by subscripts £j. 

^We do not consider the case A^'^ = A'^ — 1 corresponding to form factors of a~ , since the last ones can be 
obtained from the form factors of cr+. 
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Definition 2.1. The state parameterized by the real solutions to (j2.1ip with ij = j, j = 
1, . . . ,N^, is called a-twisted ground state in the N,^ sector. 

In the thermodynamic hmit, excitations corresponding to real set of solutions to the twisted 
Bethe equations above this a-twisted ground state can be described in terms of particles and 
holes. Note that the states associated to complex roots of the Bethe equations (like string 
solutions) are not considered in the present article. They correspond in particular to form 
factors of bound states O El [71 HH] that will be dealt with in a separate publication. All 
this means that one deals with solutions to the twisted Bethe equations where most of the £j 
coincide with their value for the ground state: ij = j except for n integers for which ij ^ j, 
with n remaining finite in the iV^ — > oo limit. To be more precise, we fix 2n distinct integers 
hi, . . . ,hn and pi, . . . ,pn such that h^ G {1, . . . , N^} and pk ^ {1, . . . , A^^^;}. The integers hk 
represent "holes" in respect to the distributions of integers for the ground state in the iV^ sector, 
whereas the integers pk represent "particles". In other words, the excited state is described by 
a set of integers {^j}^;^ such that £j = j ioi j ^ hi, ... ,hn, and ih^ = Pk- 

The terminology of particles and holes takes a clearer interpretation in terms of the counting 
function associated with the A'^^-excited state: 

^^(w) = —p^{uj) \- y ^{uj -fie,) + ^^^^ - ^ . (2.12) 

k=l 

One can argue that the counting function is monotonously increasing on the real axis. Thus 
it is possible to define unambiguously a set of real parameters as the unique solutions to 
(/ij) = j /M. One can identify, among this set of parameters, the solutions of the twisted 
Bethe equations with integers £j as being as it should be. In this picture, an excited state 
no-longer corresponds to the set of A^^ consecutive solutions fij to the equation {^j) = j /M, 
but is rather obtained by removing the solutions (m^^) = ha/M and replacing them by the 
solutions ilJ-pa) = Pa/M. It is in this respect that f^h^ stands for the rapidities of the holes 
and that fip^ stands for those of the particles. From now on, we agree upon Hj being the 
solution to ifJ-j) = j/M with .^^ being given by (|2.12p . 

Similarly as for the ground state in the A^ sector, one can also define the discrete density 
of the excited state by Pk{uj) = It is easy to see that the thermodynamic limits of this 

new counting function and its associated density coincide with the limits defined above in (|2.5p , 

An excited state is most conveniently characterized in terms of the shift function F{oj) = 
F(a;|{/ip}|{/i/i}) defined by 

F{uj) = M{i{u)-U^)). (2.13) 

The shift function describes the spacing between the root \j for the ground state in the A^ 
sector and the parameters pj defined by (fij) = j/M: 

where F{X) is the thermodynamic limit of the shift function. Recall that we consider the 
excited states in the A^^ sector with A^^ = A^ (for form factors of a^) and N,^ = N + 1 (for 
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form factors of cr"'"). Respectively we should distinguish between two shift functions F^^\X) 
and -F^'''^(A) corresponding to these two cases. Generically the shift function -F(A) satisfies the 
integral equation 

' r n " 

F{\)+ I K{\-^^)F{p)'^ = a+^- i-t^^-A +i_^[^(A-^pJ-^(A-/z,J]. (2.15) 

where 5N = N - N,,. Using ^Ol) and (f2?8]l we conclude that 

f(^)(A) = aZ(A) + ^</)(A,;UpJ - J^</>(A,^,J . (2.16) 

k=l k=l 

and 

F(+)(A) = f a - - Z(A) + ^</.(A,MpJ - 5^</.(A,/z,J + 0(A,(?) . (2.17) 
^ ^ fc=i fc=i 

3 The main results 

We now formulate a precise statement of the results obtained in this article concerning the 
thermodynamic behavior of form factors. In fact, for later use, we will give this thermodynamic 
behavior directly for the products of form factors appearing in the spectral expansion of the 
correlation function (II. 3p . However, it will be made clear that our method allows also to obtain 
the behavior of each individual form factor. The proof of the results we present below will be 
given in the next sections. 

We consider the form factors (jl.2p of the local spin operators a^, s = +,— ,z, between 
the ground state | ipg ) and some (different) excited state | V'' ) of the HamiltoniaiJl]. Here 
and in the following, we set | V'g ) = | V'({^}) )) {''^i) • • ■ > ^n} being the solution of the Bethe 
equations (j2.2p describing the ground state of the Hamiltonian, and \ = limK^-i | V'k({a*}) ) 
for , • • • , //£jv^ } being a solution of the K-twisted Bethe equations (j2.1ip with n particles 
and n holes. We will be specially interested in the limiting case where, in the thermodynamic 
limit, all rapidities {Hpj}^=i and {^^hj}^=l of particles and holes condensate on the Fermi 
boundaries. Such excited states have a zero excitation energy in the thermodynamic limit; we 
expect that they will produce the main contribution to the asymptotic behavior of correlation 
functions (this fact is already apparent in the predictions based on the Luttinger liquid theory 
and Conformal Field Theory for the long-distance asymptotic behavior of two-point functions 
[371 138| [39l Sni SH SSI HSI E]). We will also consider the opposite limiting case, namely when 
all particle/hole rapidities remain at finite distance from the Fermi boundaries. 

In the thermodynamic limit, up to uniformly ■ logM corrections, the products of such 
form factors behave as 



(3.1) 



a=0 



^The case when li)' ) = \ i>g} is trivial: T^^J^^ (m) = 0, J^^^^^^ (m) = 



2D ~ 1. 
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Let us explain these formulas. These products decrease as some negative power of the 

size M of the system, the exponents O^z, 0-\ being specified below. Due to the translation 

invariance of the model, all the dependence on the lattice spacing m — m' is absorbed into a 
phase factor, with 

N n 

j=l j=l j=l 

being the momentum of the excited state relative to the ground state. We have split the 
constant in front of the M and (m — m') dependence in two parts, Sgs' and Pg^', that we call 
smooth and discrete parts respectively. The reason for this denomination is that the smooth 
part Sss' depends continuously on the rapidities [ip. and [ihj of the particles and holes, whereas 
the discrete part P^s' also depends on the set of integers appearing in the logarithmic Bethe 
Ansatz equations ()2.1ip for the excited state. Such a discrete structure can be fairly well 
approximated by a smooth function of the rapidities [i^. and [i^. as long as they are located at 
a finite distance from the Fermi boundaries. However, as soon as the rapidities of the particles 
or holes approach the Fermi surface, the discrete structure of the form factors can no longer 
be neglected: a microscopic (of order 1/M) deviation of a particle (or hole) rapidity leads to a 
macroscopic change in T>ss'- 



3.1 Smooth parts 

The smooth parts Szz and 5^ are decomposed into 

Szz = ^ sin^ (^) . A^^ . e^^^' , = . e^^^^ . (3.3) 

Here ci'\ A'^ are functionals acting on the shift function F^^^ (A) given in SZm . while ci^\ 
aI^^ are functionals acting on the shift function F^~^\X) given in (|2.17|) . This action is such 
that the result depends smoothly on the particle/hole rapidities {/Up} and {nh} of the state 

IV'k(M)). 

Before giving explicit expressions for the quantities in (13. 3p we introduce the ivr-periodic 
Cauchy transform of the shift functions F^^^^\X) on [—q,q], 

Q 

f(^/+) {w) = — f (A) coth(A - w) dX. (3.4) 

2tti J 

-Q 

Then the coefficient Cn^ = Cn\{fip}, {fJ-h}) has the form 

n 

Ci^)[F^^)] = Co[F(^^] + 27T^Y,{F^'H^^h,-^0 + F^'H^^h,+^C)-F^'H^^P,-^C) 

-A sinh(/i/, . - fip^ - iC) sinh(/ip^ - fih, - <) ,„ 
^ smh(/ip . - fip^ - iQ smh(/ife , - fit^ - <) 
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where the functional Co[-F] reads 

Co[F] 



F{X)F{f,) 



smh (A — /i — iQ 



dX dfi. 



(3.6) 



The coefficient ci^"* = C^' {{fip}, {fih}) has similar representation, but it contains several 
additional terms 



sinh(^/j . -q-iC) sinh((U/, - g + <) 



+ ^log 



sinh(/ip . - g - iC) sinh(;Up . - g + 



(3.7) 



The functional An^ = An'{{fJ.p}, {fJ-h}) has the form 



sm vro; 



sin7rF(^)(— (7) 



2^iF(-){-g+iC) -Q sinh(g + fih^ + O 



^^^^ sinh(g + /ip^ + <) 

det [/+ 2i.^7(^)(^/;,^^;')] 



det [/ + ^K] 



(3.8) 



The equation (jS.Sp contains a ratio of Fredholm determinants. The integral operator I + -^K 
acts on the interval [—q, q], whereas the integral operators I~^^U^^^ acts on a counterclockwise 
oriented contour Tq surrounding [—9,(7]. This contour is such that it contains all the ground 
state roots and no other singularity of the kernel. The integral kernel is defined by 

U^^\w,w') = -c^(^)e^"(^'^'(-)-^'^'(-+<)) ^-^^ ~ " ^-j"^ ~ , (3.9) 

with 

/ X 1 / .^x 1 / .^x X TT sinh(tt; — u„, ) sinh(tt; — u/j, + zC) 

i^«(u;) = coth(u; + <) - Kcoth(u; - iQ, ^w) = TT . . , ^ " , . 

smh(w - /i/j^) smh('u; - + zC) 

(3.10) 

The coefficient = A^\{fip}, {/i/i}) also is proportional to a ratio of Fredholm deter- 
minants 



sin(^ 
2ttk 



e-2-P^^^(T) " sinh(/.,, - f ) det [/ + ^U^+){w,w')\ 



n 



(3.11) 



sinh(<7 - f ) J-i^ sinh(/.p, - I) det [/ + ^K] 

Similarly to the integral operator defined in (|3.9p . the operator / + ■^^U^'^^ acts on the contour 
Tq surrounding [—q,q\^ with a kernel 

[F] (wM) = Hw) ^^-^M^-g) e2-KFWW-FW(.+.c)) ^-(^'^;) , (3.12) 
^ ' ^ ^ ^ sinh(t(; - g + 1 _ g27riF(+)(«;) ' ^ ^ 



where 



sinh(?i;' + ^] 



K;sinh(tt;' — ^) 



sinh(ti;' — y) sinh(tt; — w' — sinh(?i;' + y) sinh(tt; — w' + iC) 



(3.13) 



and ^{w) is given by (|3.10|) . 
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3.2 Discrete parts and exponents 



The values of the discrete parts Vzz and , as well as of the exponents Ozz and 0^ , depend 

on whether particles and holes are on the Fermi surface or separated from it. We specify their 
expressions in the two particular cases we mentioned above: when all particle/hole rapidities 
remain at finite distance from the Fermi boundaries, or when all particle/hole rapidities collapse 
on the Fermi boundaries. Other intermediate configurations can be obtained along similar 
techniques. 



3.2.1 Particles and holes away from the Fermi boundaries 

In the first case, the exponents in (|3.ip giving the algebraic decay of the form factor with the 
system size read: 



3zz = 2n+[FX') +[F_ 

2 



.{+) 



(3.14) 



e+- = 2n+ [F'X' + 1 j + ( 

where F^f^^^ = F^'l+\±q) are given in terms of the shift functions defined in (j2.16p and (j2.17p . 
The discrete part Vzz is given by 



Vzz{{^lp],{^^h}) 



det 



sinh(//p^. - ^lhk) 
q 



- 2 n 



sm' 



J] exp i 2 /f(-)(A) [coth(A - /Up J - coth(A - fin,)] dX \ . (3.15) 



k=l 



-q 



Here a functional P*-^^ [F] has the following form: 

G2(i _ F_) G2(1 + F+) (27r)^--^+ 



[p{q) smh{2q)](F-r+iF+r 
where G{z) is the Barnes function satisfying G{z + 1) = T{z)G{z), and we have set 
q 



Ci[F] 



F'i X)Fif,)-FiX)F'{j,) 
2 tanh(A — fi) 



dX dfi 



-q 



F+ - FjX) 
tanh((; — A) 



dA + F_ 



F_ - F(A) 
tanh((; + A) 



The discrete part P-i has the form 



1^+-{{f^p}Af^h}) 
n 



det 



1 



sinh(/ip, - fihk) 



sm' 



'M+)(/.,J) 



\ /)(^;, J /3(/i, 



'PkJ 



k=l 



"sinh(^p^ - q)' 


" J 


w 


exp < 


k -q 





(3.16) 



dX . (3.17) 



exp I 2^F(+)(A)[coth(A - /UpJ - coth(A - fih,)] dX \ , (3.18) 
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where the functional reads 



IF] , XX') IF] . r"'f "'r'^u.^j;' ''''P U f \, "^l 4 p.") 

^ ^ ^ ^ sinh(2g)]i+2^+ ^| 7 tanh(g - A) ( ^ ' 



3.2.2 Particles and holes on the Fermi boundaries 



Definition 3.1. Let an excited state contain n particles and n holes. It is called critical excited 
state, if all the rapidities /ip, fih = the thermodynamic limit. We also say that this excited 

particles, resp. 



state belongs to the Pf class if it contains particles, resp. holes, with rapidities equal to 



±g such that 

Up — = n'f^ — Up = r, r G Z. (3.20) 
The corresponding form factors are called critical form factors of the P^ class. 

The rapidities of the particles and holes in an excited state belonging to the Pj. class are 
all located in a close neighborhood of itg. As a consequence, it is useful to re-parameterize the 
integers describing the position of particles and holes according to 



Pj =pj + if Ptpj =q, 

(3.21) 



= 1 - Pj , if Ptpj 



hj = Nf, + 1- if /x/j^ =q, 
hj = hj, if fj^h. = -q. 

All the integers {p^} and {h^} introduced above are positive and vary in a range such that 



lim = lim = , (3.22) 

which means that fip and fih indeed collapse to the Fermi boundary in the thermodynamic 
limit. 

For and form factors belonging to the P^ class, the discrete parts are given by 



G^(l + Ff )G^(1-f1-^) / sin(7ri^,g) y"'- j^ sin(vrFg) ^"'^ 



vr / \ vr 



and 

,G2(2 + F|+^)G2(1 -F1+^) /sin(^F;+0\"'"' /sin(^F, 



X Rnt,nti{P^hmK>) K~,^-i{p-},{h-}\-F'y) , (3.23) 



G2(2 + F^+^)G2(l -FJ+)) V TT J y TT 

Rni,ntiiP^hm\l + Fi^^) K^,^~i{p-},{h-}\ - fW) . (3.24) 



Here we agree upon F^^^^X) = F(^/+)(A) + r and -F^fi"^'' = Fi:''^'^\±q). The functionals 
p(^) and are given respectively in (|3.16p and (|3.19p . The coefficient Rn,m{{p}, {h}\F) is 
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defined as 



n m 
^2 



uipj-pkruihj-hkr . 

Rn,m[{p}Ah}\F) = r r 1 rz, 1 I' (3-25) 

n UiPj+hk-ir V 

j=i k=i 



{Pk} , {hk} 



where we have used the standard hypergeometric type notation for ratios of T functions: 

\ i j 

' =l[T{ak)-l[Tibk)-\ (3.26) 




k=l k=l 



The exponents 9 in this case are 

(3.27) 



l+FW)^(FW)^ 

Note that the shift functions i<'(^/+)(A) enter the exponents ()3.27p only in the combination 
fI:''^^\x) = F(^/+)(A) + r. It is easy to see that fJ^''+^(A) are nothing else but the shift 
functions of the (a + r)-twisted ground states (see Definition 12. 1|) in the N and + 1 sectors 
respectively. 

Remark 3.1. One can write down the results of this section in a more symmetric form. For 
this we introduce two constants /zp describing the shift of the utmost parameters /ii and 
of the (a + r)-twisted ground state with respect to the utmost parameters Ai and Aat of the 
ground state: 



/_= hm Mp(Ai)(/ii-Ai), 

N,M^oo 
N.M^oo 



(3.28) 



(z) -L. 

Then it is easy to see that f± = F.^J. for the form factors. However for the a'*~ form factors 
(where Nf^ = + 1) one has /_ = -FJ^\ but /+ = 1 + F^~^ . Then in both cases the exponents 
([3:271) can be written as 6 = fl + fl. 

It is crucial to compare the results (|3.27p with the critical exponents appearing in the 
asymptotic behavior of the two-point correlation functions predicted in [371 EH EH SQl El H2l 

ISIIMIIMIIMIES]. 

In the case of (9^^(r), it follows from (fTTUD that ^^''^(A) = (a + r)Z(A). Therefore, 

for a = 0, we obtain that 9zz{f) = 2r^Z^((7), \r\ = 1,2, ... . These numbers coincide with the 
critical exponents appearing in the asymptotic behavior of the (cfcr^^^) correlation function 
and are associated to the oscillating term with momentum 2rkF- 

The exponent 9^ (r) has similar properties. In this case, F^^\X) = (a + r- 1/2)Z(A) + 

(/)(A,g). It follows from equations <\2M and (I2.10p that 

Z{q) + Z-\q) 



1 + (t){q,q) 



2 

Z{q)-Z-\q) 



(3-29) 
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From this we find 

Fr,+ + l = {a + r)Z{q) + ^^, 

y-x, . (3-30) 
Fr,. = {a + r)Z{q)-^^. 

Thus, at a = we have 9^ (r) = Z~'^{q)/2 + 2r^Z^(g), \r\ = 0, 1, ... . Again, these numbers 

coincide with the critical exponents appearing in the asymptotic behavior of the {<Tfa^_^_^) 
correlation function and associated to the oscillating term with momentum Irkp. Thus, we see 
that the behavior of critical form factors with respect to the size of the system coincides with 
the asymptotic behavior of two-point correlation functions with respect to the lattice distance. 



4 Form factors and scalar products 

The solution of the quantum inverse scattering problem enables us to express the local spin 
operators in terms of the entries of the monodromy matrix [10\ |23] : 



< = T^-\-Km ■ tr (T(-</2) a') ■ r-™(-iC/2). (4.1) 

In the l.h.s. of this expression, the symbol (s = it, z) denotes the corresponding local spin 
operator at site m, whereas the symbol a* appearing in the r.h.s. should be understood as a 
2x2 Pauli matrix multiplying the 2x2 monodromy matrix (12. ID . 

Using ()4.ip . one can reduce the computation of the form factors of local spin operators in 
the finite XXZ chain to the one of the scalar products [10]. We have 

JV JV 

/.^ i(m-l) J2 [Po{lJ.ej)-po{^j)] / i E [po(m«,- )-Po(Aj)] 



da7ri\\i;^m)\\-\mX})\\ 



, (4.2) 

Q! = 



(z) 

which is nonzero only if N,^ = N. The form-factor J-"^ ^, is obtained by complex conjugation 
^Ipg^iim) = (^-7^^^^g("T-)) • Similarly, one obtains 

JVk JV _ 

J-(-) im) = f-n^+-+i ^^("^-i)EpoK)-^^Epo(a.) (^.(M)|i?(-f)|V^({A})) 

^'^^ a{-^i)\\M{m■\mm\ -o' 

JV JVk 

"^'^ a(-4)llV'.(M)MlV'({A})|| -o' 

which are nonzero only \i N^ = N + 1. In (Oj) and (fOl) . a{v) = sinh^^ {v - iC/2) denotes the 
eigenvalue of the operator A{y) on the reference state which, in the case of the spin-1/2 chain, 
is the ferromagnetic state with all spins up. In the following, we will also use the notation 
d{i') = sinh*^ {y + iC,/2) for the eigenvalue of D{u) on this reference state. 

Remark 4.1. Although it is possible to set a = directly in (j4.3p and (j4.4p . we prefer to take 
the limit a = only in the very end of the calculations to study the properties of the scalar 
products for general a. 
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In all the above scalar products, one of the states is an eigenstate of the twisted transfer 
matrix. There exists an explicit determinant representations for such scalar products and the 
associated norms: 

Proposition 4.1. Let {X]i satisfy the ground state Bethe equations (|2.2p and {/u}^" solve 
the system of a-twisted Bethe Ansatz equations (j2.1ip . Then the following representations for 
the scalar products hold: 



-LJ- smh(Aa - nej -^-^ ' ' 



a,b=l 
1 - K 



„2i7rF(Aj) 



N 



^ sinh(g + A.-zC) ^^^^ 



ZTTl 



, (4.5) 



N 

B{-iC/2) m{X})) = 5n.,n+i ■ ai-iC/2) sinh (-<) H " e^'^^^^^^] } 

N 

n sinh(Ab - iC/2) n+i ( N . , . ^ 

b=i TT j J/ N TT sinh(/i£^ - h - <) 



Af+l 

n 

6=1 



n sinh(/X4+</2) -1 >^ b=i "-^'^ ^> 



(4.6) 



/n these expressions F{X) (j2.13p depends on whether we consider the N sector or the {N + 1) 
sector (see ()2.12p ). The integral operators I + ^U^^^ and I + ^U^'^'^ act on a closed anti- 
clockwise oriented contour Tq surrounding the interval [—q, q] where the ground state roots 
{Xa}a=i condensate and containing no other singularity of the kernels. The last ones have the 
form 

sinh(u; - Aa)sinh(u; +<) i _ ^2TTiF{w) ' 

^ ' ^ ^J-^ smh{w - fie^ + %() ^J^ smh{w - Xa) 1 - e^^^^i^) ' 
where and IC^ are given respectively by (|3.10p and (j3.13p . 

The representation (14. 5p was obtained in [JU [25l [ini [36] . The representation (14. 6p is derived 
in appendix iBl 

We also recall the finite size determinant representations for the norms of Bethe states: 
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Proposition 4.2. JFJ[ EE IS^ Let ne^,. . . satisfy the system (j2.12p . Then 



n sinh(/x£„ - fie,) 



a, 6=1 



xdetej.^^ (4.9) 



One can see from another representation for the scalar product (|B.ip that, when SJq = 



(V'.(M)IV'({A})) = (V.(M)I^({A})) 



(4.10) 



(4.11) 



and 



(V^.(M)|i?(-|)|V({A})) 



_ (4.12) 

For the calculation of the two-point correlation functions (cr^cr^i) (resp. (c^cj",)), we 
need actually to sum up the products -^^^^^/("i) -^^^^^("2.') (resp. -^^^^/("i) J-"^;J^(m')) over all 
eigenstates \ip' )■ We have the following result: 

Proposition 4.3. The products of two form factors can be written in terms of the former 
scalar products as 



2 sin 



2 r,. 



S 



N 



a=0 



a=0 



where 



S 



N 



s 



N 



(V^.(M)|i?(-f)|^({A})) 



a2(-iC/2) 



^'.(M)MIV'({A})|| 



(4.13) 
(4.14) 

(4.15) 
(4.16) 



In (j4.13p and (|4.14p . Vex denotes the relative excitation momentum of the excited states in the 
sector in respect to the ground state in the N sector, 



Vex = ^PoilJ-i,) - ^Po(Aj). 

i=i i=i 



(4.17) 



Remark 4.2. Here, we did not consider the product J^^ ^'i^) -^^^ ("t-Oj since the equal-time 
correlation function {(7^^^') '^^^ t)e obtained from (i7+(T~,) by the replacement m ^ m' . 
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The calculation of the thermodynamic limit of Vex is a relatively simple problem. Indeed, 
it follows from (plil) . (l2l^ that 

3=1 j=l k=l 

(4.18) 

Here we have used the antisymmetry of the bare phase A) = and S,K,ij^e ) = ^j/^-, 

£,{Xj) = j/M. Then, using p.6p . we obtain that the thermodynamic limit Vex of Vex is given 
by the expression (13. 2p . 

The non-trivial problem is the computation of the thermodynamic limit of the scalar prod- 
ucts Sf^ and S^. This problem will be treated in the remaining part of the article. 



5 Thermodynamic limit of 

In the paper [I] we studied the thermodynamic limit of the scalar product Sf^ in the particular 
case where the state IV'kIIa*})) was the a-twisted ground state in the sector. In the general 
case, when the excited state contains n particles and n holes, the result also depends on their 
rapidities {/ip} and {f^h}- As we have already announced, such dependence is not sufficient 
to characterize completely the thermodynamic limit of the form factor. Indeed, the latter 
decomposes into a product of a smooth and a discrete part. The smooth part can actually 
be completely described in terms of the rapidities of the particles and holes. However, for the 
description of the discrete part one should use also the integer numbers {pa} and {ha}- 

In the first part of this section, we explicitly factorize Sf^ into the aforementioned prod- 
uct. Then, we investigate the thermodynamic limit of the smooth part, postponing the more 
complicated analysis of the discrete part until subsection 15.31 

Note that in this section we deal only with the excited states in the N sector. Therefore the 
thermodynamic limit of the shift function -F(A) is equal to F^^\\) given by (|2.16p . However, 
in order to lighten notations we omit the superscript (z) throughout this section, denoting the 
limiting value of the shift function simply by F{X). 



5.1 Representation of the scalar product S 



N 



Assume that the excited state |V'k({a*})) is an excited state in the sector containing n particles 
with rapidities {fip^}"^^-^ and n holes with rapidities {fJ'ha}a=i- ^ne can fairly expect that the 
limiting value of Sf^ depends on the thermodynamic limit of these rapidities {fj.p} and {nh}, 

i.e. limSf^ = S^{{Hp},{Hh})- 

Using Propositions 14.11 and 14.21 as well as ()4.15p , it is readily seen that has the following 
representation: 

SUM, M) = A^n\M^ M) ■ DPi{f,p}, {/.,}) ■ exp {ci^Hifip}, M)} , (5.1) 
where 

c^\M, {«}) = E log ""l;!>°"?"n "w""'""'"'.? ' P-^' 

smh(Aa - Xb- <) smh{fie^ - ^4 - iQ 
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D'^^\{fip},{fih})= [det 



2 N 



n 



sin^ 7rF(Aj 



(5.3) 



and 



N 



sinvra -r-J smh.{q + Xa — iC) 



n 



detr„ 



detN eg^ -detivej^ 



(5.4) 



simrFi-q) smh{q - - iC) 
We remind here that U^^"^ has been defined in (j4.7|) and that the N x N matrices 0^-^'' and 0^-^^ 

(5.5) 



are 



5.2 Thermodynamic limit of the smooth part 



27rMp«(^4) • 



The smooth part of the scalar product consists of factors A^^^ and Their hmits can be 

computed exactly in the same manner as in [Ij. We set 



d'\{pp},{f,h}) = ^ lim c5;)({Mp},{^4), 

A',i\/— J-oo 



(5.6) 



(z) (z) 

We illustrate the main idea behind the calculation of Cn and An on the following toy- 
example. Assume that /(A) G C^(M), and consider the computation of the thermodynamic 
limit of the sum 



N 



To compute this limit one can decompose the sum as follows 



(5.7) 



N n n J. 

Sf = Y.[f{^^J)-f{XJ)]+Y,[fi^^P,)-f{^^h,)]^Y.ifM-f(^'h,)]+ F{x)f'{X)dX. (5.8) 

i=i i=i i=i -q 

There we have used the leading behavior in M^^ ()2.14p for the spacing between pj and Xj. It 
is convenient to introduce a modified shift function by 



i^w(A) = F(A)X[_,„j(A)+5; 



(5.9) 



where X[a,b] is the characteristic function of the interval [a, b]. Then equation ()5.8p can be recast 
simply as 



Sf^j F„^od(A)/'(A)dA. 



(5.10) 



Therefore, provided that one extends the integration contours to the whole real axis and replaces 
F with Fmod^ one formally reduces the computations to the case considered in [I] where one 
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was dealing with an excited state IV'kCI^}) having no particles or holes. From this observation, 
the results of section [3] follow straightforwardly (by a mere replacement of the shift function by 
Fmod) from those of [Ij. In this way, starting from the equations (I5.2p . (15. 4p we arrive at ()3.5p . 

das]). 

Thus, the limits An and CA are well defined for arbitrary positions of particles and holes. 
They depend on {/^p} and {fih} only. In particular, they do not depend on the underlying 
integer numbers. 

At this stage it is interesting to consider the limiting case where, in the thermodynamic 
limit, all rapidities of particles and holes condensate on the Fermi boundaries. 

Corollary 5.1. Let |V'k({/^})) belongs to the class as in (|3.20p . and denote 
eg = Ct\{+q}^. U {-q}^.,{+q}^,^ U 

where subscripts show the number of elements in the corresponding subsets. Let Fr{X) = F{X)+r 
and Fr{w) be the iir-periodic Cauchy transform of Fj.{X). Then the coefficient Cn). takes the 
form 

Cg = C7o[F,], (5.12) 
where the functional Cq [F] is defined in (13. 6p . The coefficient An}- becomes 



smvra 
simTF{—q) 



2 



g-27riF,(-g+iC) 



det [l + ^^U^'\w,w')] 



det [I + ^K] 



(5.13) 



where 



Jjiz). ^ ^2.i(Fr(w)-FAw+i()) K^jw - w') - K^{-q - w') 



(5.14) 

Proof. This follows straightforwardly from Fmod{^) = Fr{\) X[-q,q\{^)- D 
Using (US]), (fZTB we find that 

Fr{\) = {r + a)Z{\), (5.15) 

and we see that up to a replacement a + r ^ a our results coincide with the ones obtained in 
[T]. This coincidence is not accidental. 

In the thermodynamic limit, the class contains an infinite number of states. Indeed, 
knowing that a particle or hole's rapidity equals to itg in this limit does not allow to fix the 
corresponding integers p (resp. h) unambiguously. For instance, knowing that fJ-p = q only 
allows one to say that p/M D. The latter is satisfied as long as one chooses p = N + um 
with Um = o (M), but arbitrary otherwise. 

It is clear that quantities having well defined thermodynamic limits (such as Cn,r and An'r)^ 
must only depend on the macroscopic realization {^p} and {/i/i} of the excited state and not 
on the microscopic quantum numbers {p} and {h} leading to such a macroscopic configuration 
of the rapidities. Therefore, the calculation of their thermodynamic limit in the case of an 
excited state belonging to the P,- class can be done by choosing any of its representative. One 
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of these representatives is given by the (a + r)-twisted ground state. Indeed, set = = 
and choose the integers describing the holes and particles according to 

hk = k, pk = N + k, k = l,...,r. (5.16) 

The excited state is thus described by the set of integers, 

(j = j for j = r + 1,. . . ,N, and £j = N + j for j = 1, . . . , r . (5.17) 

Observe now that if we replace a by a + r in the r.h.s. of (12. lip , then up to a re-ordering of 
the equations, we obtain the same set of integers (j5.17p . Thus, the above representative of the 
class Pr coincides with the (a + r)-twisted ground state. Therefore the smooth part of all the 
critical form factors belonging to the class coincides with the one of the form factor of the 
operator taken between the ground state and the (a + r)-twisted ground state. 



5.3 Thermodynamic limit of the discrete part 

As we have seen, the modified shift function (|5.9p can be used for the calculation of the ther- 
modynamic limits of sums (products, determinants) in the case when we deal with smooth 
functions of the rapidities of the excited state. One cannot apply this method in the case of 
the factor D^^\{np} , {Xh}) as it depends on the Cauchy determinant of {A} and {^}. In the 
thermodynamic limit, certain entries of this matrix become divergent and more care is needed 
for the calculation of this limit. 

(z) 

Proposition 5.1. The discrete part Dj^ behaves in the thermodynamic as 



k=l 



The factor Eq depends explicitly on the integers {ha} and {pa} which are the quantum numbers 
defining the excited state 

n 

n {hj - hk){pj - Pk) vilJ-h, , IJ-hk) ^ilJ-p, , IJ-Pk) 
Eo = n > (5-19) 

n (pj - hk^^'^if^pj^i^hk) 
j,k=i 

where ip (A,/i) = (27rsinh(A — /x)) / (p(A) — pin)) and p{X) is the dressed momentum ()2.6p . The 
factors Pn and Hn also depend on the quantum numbers defining the holes and the particles 

^ ' ^ e^[^'^']K)vrV(/^hJ V hk,N + l-hk J ^ ' 
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There we have introduced the functional J which is defined in terms of the dressed momentum 
p{X) and acts on the shift function F: 

q 

J[F]H = 2 1 |f (A) dx log y(A, ^) + ^ ~ ^^^^ p' (A) I dX . (5.22) 
Finally, we remind that the functional T>^^^ is given in (|3.16p 

Remark 5.1. The definition of Pm {iJ-p^.^Pk) is given in the case where the rapidities of all 
particles are to the right of the Fermi zone { ie pk> N). If some particles have their rapidities 
to the left of the Fermi zone {pk < 0), then the corresponding arguments of the P- functions in 
(I5.20p become negative integers. The aforementioned formula remains however valid provided 
that the arguments of these F-function are understood as limits 

^^''^ -lim,,i>^ = (-l)-i:i^±l^. (5.23) 



T{pk-N) e^oT{pk-N + e) ' ' r(l-pfc) 



Proof. Following the strategy applied in [T], we multiply and divide the original Cauchy deter- 
minant by the Cauchy determinant of the counting functions ^k(A). Let 

_ sinh(A-/i) _ sinh(A-/Li) /r 

UA)-^(^) C«(A)-^k(//) 

(z) 

After some algebra, we recast D)^ {{fip}, {Xh}) into the following product: 

n n 

D^NiM^ = Eo llHk llPk Dn,o . (5.25) 

k=l k=l 

Here Eq depends only on the particle/hole rapidities and the corresponding integers: 

n 

n (hj - hk) (pj -Pk)^K ifJ-hj ifJ'hu)'?^ ifJ-p, > Ptpfc ) 
Eo = — E • (5-26) 

n (pj - hkY(p'^{iip^,iihk) 
j,k=i 

The factors Pk depend on the rapidities ^p^, integers pk, the ground state parameters {A}, and 
the roots //j such that ^^(/^i) = j/-^; j = ■ ■ ■ , Their explicit representations are 



N 



,2 



Ml^Pk) f}^ \ j-Pk- F{Xj) ^« (PPk^^j) J \Pk + F{fip^) , pk-N J ' 
The representations for Hk are similar to (j5.27p . but they depend on and hk'. 
^ sin^(7rF(/xfeJ) f j - hk - FjXj) _ (p^ {l^nAj) 1 



[ hk, N + l-hk ■ y ■ ' 
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Finally, Dj\[fi does not depend on the particle/hole rapidities and the corresponding integers 



(5.29) 



The fact that Eq — )• £'0 follows from (p^ — )• ip. The calculation of the limit of the first 
product in (15.290 is based on the definition of the shift function (12.130 : 



-n(l-|i;)) -exp(-/nA).Aj ^e--.. (5.30) 

The limit of the remaining part of D^^q was computed in [Ij, what gives us 

Dn^o = M-n-^'- . P(^) [F] . (1 + O (i^^) ) . (5.31) 

Finally, the limits — > PN{fip^jPk) and Hj^ — > H^^^jLi^^^hj^) follow from 

N ( ^ ^ ^ 

hm TT l^i^ZM|4^ =e^[^l(..). (5.32) 
iv,M^oo 11 [j -k- F{Xj) 'fK (Mfc, ) J 

This formula is proved in appendix |A] (see (lA.SO ). □ 

Assuming general positions of the particles and holes, we can not make further simplifica- 
tions in the formula for -D^^- Therefore we now study two limiting cases of interest, exactly as 
we did for the smooth parts: first when all particles and holes are separated from the Fermi 
boundaries, and then when all particles and holes are in the vicinity of the Fermi zone (critical 
form factor in the class P^)- 

Corollary 5.2. Suppose that, in the thermodynamic limit, particles and holes are separated 
from the Fermi boundaries: fip^ 7^ ±q and fi^a ^Q- Then the thermodynamic limit of the 

(z) 

factor becomes a smooth function of the rapidities {fJ,p} and {fJ-h}' 

= M---^^-n . Di:\M, M) . (1 + O (i^^)) , (5.33) 

where Dn\{fip}, {fi^}) is given by the expression (j3.15p forV^z- 
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Proof. If the rapidities of particles and holes are separated from the Fermi boundaries, then all 
the arguments of the F-functions in ()5.20p and ()5.2ip are large and we can apply the Stirling 
formula for their simplification. Then 

lim F^ I ' Pt'~ ^ ^if^Pk) ' + 1 - /ifc - F{fihJ > + F{^hJ \ 
7V,M™(^ \ Pk-N , pk + F{npJ , N + l-hk , hk J 

lim ' ' 



( pk-N V 
\ Pk ) 



N,M^oo \ Pk J \ hk 

Pi^^Pk)-PiQ) V^^^'"^ f PiQ)-pM V'^^'^'^^^ 



MfJ'Pk) - Pi-Q) J Kpii^h^,) - p{-q) J ' 

where we have used (j2.6p . Combining (j5.34p with the explicit form (|5.22p of J (u;), we arrive at 

sin2(7rF(^/,J) 



Pn {^J'Pk,Pk) Hn ifj-hk^hk) 



■^'^pifj'hJpifj'Pk) 

X exp I 2 /F(A)(coth(A - /Up J - coth(A - fih,)) dX \ . (5.35) 



The limit of Eq is almost trivial: 

lim ^oM2" = f det — — ^ . (5.36) 

N,M^oo \ n Smh(/ip^ -^/ifc)/ 

Gathering all these results, we get the claim. □ 
Hence, in the limit when the rapidities of particles and holes remain at finite distance from 

(z) 

the Fermi boundaries, the discrete part Dj^ , up to an M-dependent normalization, can be 
approximated by a smooth function of these rapidities. However, as soon as some rapidities 
go to ±g, then certain integrals in the second line of (I3.15P become ill-defined. There can also 
appear singularities in the Cauchy determinants present in the first line of ()3.15p . Therefore, 

(z) 

the discrete structure of manifests itself when particles and holes condensate on the Fermi 
boundaries. 

We now describe the asymptotic behavior of an excited state \ipK{{ni^}) belonging to the 
P,. class. This means that the rapidities of all particles and holes are located in vicinities of 
the Fermi boundaries with the condition (|3.20p . Any such state, can be parameterized by a set 
of integers as given in (I3.2ip . 

(z) 

Corollary 5.3. The asymptotic behavior of Dj^ for an excited state in the class parame- 
terized as in ()3.2ip reads 

d'^\{+<iU U {-l}n,d+<l}nt U {-cz}.-) = M-^^.-^^ (l + O (i^^)) , (5.37) 

where DnJ- is given by the expression ()3.23p for Dzz ■ 
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Remark 5.2. Thence, when particles and holes are in the vicinity of the Fermi boundaries, 
the behavior of DVr cannot be called a thermodynamic limit. Indeed, the factors R+ + and 

jv rip 



-R^- thus the form factor, depend explicitly on the microscopic characteristics and 

rip 

of the excited state. 



Proof. In this case, the limits Eq, PN{np^,pk) and Hi\/{fj,h^, hk) should be re-calculated. 

Let us first consider the coefficient £"0- Using that f{q,q) = ^{—q,—q) = p~^{q) and 
"Pil, -q) = ^{-q,q) = sinh(2g)/L», we get 

^ / D / 1 \^ 
lim Eo = lim det — . (5.38) 

N,M^oo \p{q) smh.{2q) J N,M^oo \ n pj - hk J 

Using the specific parametrization (I3.2ip of the integers belonging to the P,- class we see that 
the thermodynamic limit of the Cauchy determinant in (I5.38P can be factorized in two parts: 

with 

n m 
\[ij>3-Pk?W{h,-h,f 

LuAip]^ m = '-^^^^hr^ — ■ (5-40) 



n uipj+hk-ir 
j=i k=i 



This leads eventually to the following estimate 

\2r2 2n 



lim £;o(Mp(g)sinh(2g))^'^ = L + L ({p"}, (5.41) 

A',Ai— >oo P ' h P ' h 



When computing the limits of Pk and H^, some of the arguments of the F-functions are not 
large anymore. Therefore, we can only use Stirling formula partly. We have 



k=l 



2 ( Pk , Pk - N + F (fipj , N + l-hk- F {nhj , hk + F (/i/,^ 



Pk- N , Pk + F {fj,p^) , N +l-hk , hk J 

^-2.(F++F_) p2 ( {P+ + F+} , {h+ - F+} , {/i- + F_} , {p- - F_} 
i , {h+} , {h-} , {p-} 



(5.42) 



We precise that the notation {p+ + means {pj' + F^}^^^ and similarly for other sets of 
parameters in (I5.42p . Combining all these results we obtain 



lim ^0 {Mp{q) sinh(2g))2'^'+2KF++F_) tt p^Q^ 



k=l 

2n 



(^) fl„;,<({pn.{^ni^.)iJ„,-,„;({p-).{^-}l-^-) 

V g -q ) 

After some simple algebra we get the claim. □ 
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To conclude this section we would like to stress that the evaluation of the Sf^ thermodynamic 
limit was done without using the explicit form of the shift function F{X) (I2.16p . We have used 
this representation only in ()5.15p in order to relate the smooth part of the critical form factors 
of Pr class with the one corresponding to the (q + r)-twisted ground state. However in all 
other respects F(X) played the role of a free functional parameter. 



6 The scalar product S]^ 

We now study the scalar product for | Tp,^{{^}) ) being an a-twisted excited state with n 
particles and n holes in the + 1 sector (i.e. N,^ = + 1). As previously, {lJ'pa}a=i 
{l^ha}a=i denote the rapidities of the particles and holes respectively. The notation F{\) now 
means the shift function in the + 1 sector, whose thermodynamic limit F^~^\X) is given by 
(j2.17p . However, like in the previous section, we omit the superscript (+) in order to lighten 
the notations. 



6.1 Representation for the scalar product S 



N 



In comparison with the Sf^ case, the excited state now depends on + 1 parameters /i. 
Using the determinant representation given in Proposition 14.11 as well as the norm formula of 
Proposition 14.21 we get the following representation for (j4.16p : 



Sj,{{f^p}, M) = A^n\M, M) ■ exp {fih})} ■ D'^N\{fip}, M). (6. 



11 



Here 



N+l 



+1) a=l 



n sinh (Aat+i - mj 



^ ' [[ smh (Aat+i - \a) 

a=l 



D'-^l,i{fip},M), (6.2) 



where -Cj^^-^ is given by ()5.3p with A^ replaced by A^ + 1. It is expressed in terms of the A^ + 1 
parameters fi£^ as well as the A^ + 1 parameters Xj defined by ^ (Aj) = j/M, ^(A) being the 
counting function (|2.4p for the ground state. In other words Aj, j = 1, . . . , A^, are Bethe roots 
for the ground state, whereas Aat+i is the point where M^(A) takes its next integer value, ie 
M^{Xn+i) = N + 1. 



The coefficient is now modified according to 



exp 



c{,+)({Mp},{^4) 



7V+1 

n 

a=l 



sinh(Aa - Aat+i - iC) 



sinh(^^^ - Atv+i - iC) 



exp 



> (6-3) 



where, up to the evident modifications stemming from A^ 
Finally the factor A^j^^ has the following form: 



N + l, Ci^l^ is given by ([El. 



sin ( 
2irK 



N 

n sinh(A« 

a=l 



2 ^ 



N+l 

n sinh(^£„ 

a=l 



2 - 



det 



w, w 



det^+1 •det^veg) 



(6.4) 
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and we refer to (j4.8|) for the definition of the kernel U^~^\w,w'). 



6.2 Thermodynamic limit of the smooth part 

The smooth part of consists of the coefficients A^^^ and expC^\ The computation of 
their thermodynamic limits and does not contain any subtleties comparing with the 
derivation described in section 15.21 The use of the modified shift function (15. 9p formally reduces 
the computations to the case where the excited state does not contain particles and holes. In 
this way, starting from the representations ()6.3p , ()6.4p we arrive at the equations ()3.7p , (j3.1ip . 
One can easily see that the coefficients ci"*"^ and A^rt^ are well defined for general positions of 
particles and holes. In particular, for the critical form factors of the P,. class the smooth part 
effectively depends on the shift function Fr{X) = F{X) + r. 

6.3 Thermodynamic hmit of the discrete part 
Proposition 6.1. The discrete part i?^^ behaves in the thermodynamic as 

n 

D^n\{^^p], {l^h}) = M-^--^^++^^" Eo {{Hh,l^p] , {ha,Pa}) ■ n ^^+1 if^h^hk) 

k=l 

n ^^^+1 if^P.,Pk) ■ VM[F] . (l + O (^)) , (6.5) 



where 



HN+iifJ'hk^h) = HN+i{f^Pk^Pk) + - hk- F{fihJ) (f (g,)UfcJ, (6.6) 
Pn+1 {f^pk,Pk) = Pn+1 {tJ'Pk,Pk) {N + l-pk- F{fip^)f(p^{q,fip^). (6.7) 

The functions entering these definitions are the same as in Proposition 15. il but one should 
replace N by N + 1 for the coefficients Hn+i and P/v+i • We also remind that the functional 
P*-^-* is given in ()3.19p . 

(z) 

Proof. As the limit of Dj\f!^i is already known, it remains to evaluate the limit of the product 

N+l 

n sinh2(AAr+i -/x^J 

P = ^ . (6.8) 

n sinh2(AAr+i - A^) 

a=l 

Harping on the steps in the analysis of D^j^\ we multiply and divide by the counting function 
^(A) (Hai). We obtain 

A^+l / ^ \2 N+l 



n (N+i-ia-Fifiej) n ^'(A^+i,M£„ 



P=— ^ ^ • (6-9) 

a=l a=l 
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Then using Corollary I A. II we find that, in the thermodynamic limit, the product P behaves as 



P 



eJlF]{q) 



N + 1-F+ 
N + 1 , -F, 



n 

k=l 



N + l-hk- F{iJhk ) f{q, f^hk 



1 2 



with J[F] {q) defined as in (I5.22p . Using now that N is large, we have 



(6.10) 



gJ[F](g)p2| 



N+1-F+ 
N+1 , -F. 



r2(l + F+)sin2 ttF^ 



exp < 2 



7r2 {p{q) sinh(2g) Mf^+ 
Taking into account (I3.16p . (13.190 we find 



F{X)-F+ 
tanh (A — q) 



(6.11) 



sin2(^F(Aiv+i)) 



k=l 



N + l-hk- F{phk) 'Pil^lJ'hk) 



(6.12) 



It remains to combine the obtained result with the known limit of D 
into (I6.2p . A few algebraic manipulations lead then to the claim. 



(^) 

N+l 



and to substitute it 
□ 

We now particularize this result to the two limiting cases we considered previously, namely 
when all particles and holes are separated from the Fermi boundaries, or when all particles and 
holes are on the Fermi boundaries. 

Corollary 6.1. Suppose that, in the thermodynamic limit, particles and holes are separated 
from the Fermi boundaries: ^ +q and ^Q- Then the thermodynamic limit of -D^"^ 

becomes a smooth function of the rapidities {/ip} and {ph}-' 



where D^\{np}, {ph}) is given by the expression ()3.18p for Vj, 



logM 
M 



(6.13) 



Proof. If the rapidities of particles and holes are separated from the Fermi boundaries, then 



n 

fc=i 



N +l-pk- Fjfip^) ^ fiq,HpJ 
N + 1 - hk - F{fihk) fiQ,f^hk) 



n 



/ sinh(g-^pj 



\ \smh{q - phk, 



(6.14) 



Substituting this limit into (16. 5p and using the results of the previous section we arrive at the 
statement. □ 

Corollary 6.2. The asymptotic behavior of D^^^ for an excited state in the P,. class parame- 
terized as in ()3.2ip reads 



D^N^ i{+lU U ^-'^K ' ^+^>< U {-q}^- ) = M-(^^'++i) ( 1 + O 

where d'^} is given by the expression (j3.24p for 



log M 
M 



(6.15) 
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Proof. In this case 



n 

k=l 



N + 1 - hk - F{nh^) ip{q,iJhJ 



1 - F^)^ 



Substituting this limit into (16. 5p and using the results of the previous section we arrive at the 
statement. □ 



Conclusion 

In this article we have provided the first steps towards the study of correlation functions via 
the form factor approach. For this purpose we have calculated the thermodynamic limits of 
the so-called particle/hole form factors. Although we have considered the specific case of the 
XXZ chain, our method is straightforwardly applicable to other massless integrable models 
solvable by means of the algebraic Bethe Ansatz and having determinant representations of 
their form factors. In particular, one can use it for the calculation of the thermodynamic limit 
of form factors in the model of one-dimensional bosons. Determinant representations for the 
form factors of this model in the finite volume were given in [W\ I2U| . Their thermodynamic 
limit can be directly obtained from the results given in the present article. 

Our results show that the idea to replace, in the thermodynamic limit, the sum over excited 
states by an integration with respect to particles and holes should be essentially modified. 
In particular, a formal replacement of the discrete sums over n-particle/hole excited states 
by integrals as 

/n 
4lV ("^) 4' i ("^') n P(^p. )^(^^. ' 

(6.17) 

leads to senseless results. First of all, even in the region where the form factors depend smoothly 
on {/Up} and {nh}i the coefficient M^" does not compensate the factors M"^*^', therefore in 
the above expression, the r.h.s. vanishes. On the other hand, the corresponding integral 
becomes divergent when the rapidities of particles and holes approach the Fermi boundaries 
(see (|3.15|) . (j3.18p ). Lastly, due to the discrete structure of the form factors when some rapidities 
agglomerate on the Fermi boundaries, one cannot replace the sum over the excited states in 
the Pr class by an integral. For such excited states, one has to take the microscopic structure 
of the excited state into account and perform the discrete sums over the parameters p^, 

In a forthcoming publication |54] . we will describe a way to overcome these difficulties. We 
will show that in the asymptotic regime (large lattice distances m between local spin opertors) 
only particles and holes having their rapidities close to the Fermi boundaries contribute to the 
form factor sums. In such a limit, one can send the rapidities of the particles and holes to itg in 
the smooth part of the form factor. Hence, it becomes a constant that only depends on the 
class 5ss'='?^s' ({i?}) {i^})- III its turn the summation over the integers {p} and {h} present 
in the discrete part of the form factors leads to the natural re-scaling M — > 27rm of the system 
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size M into the distance m between the operators. This mechanism explains the appearance 
of the critical exponents Ogs' in the asymptotic behavior of the correlation functions. 
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A Summation identities 



Lemma A.l. Let f G (^^([0, a]) for some a> D. Let (/im) be a sequence of integers such that 
^ tends to some finite value Xh £ [0, a] when M — )■ +00. Then, in the limit N,M ^ +00, 
N/M D, the following sums vanish: 

N 



lim >^ 

N.M^Qo ^ 

k=l 



n>2. 



(A.l) 



Proof. Let n = 2. Then 



N 



I/(m 



fe=i 

k^hM 



M> 



/(l 



M ' 

M , 



k—h]\ 
M 



■/'( 



M • 

M , 



M2 



fe=l 
k^hM 



V M I 



+ 



M . 



M 



N 

Y. 

k=l 
k^hM 



1 



k — Hm 



(A.2) 

The second term vanishes in the limit, as it is of order ^°|^^ . The first term vanishes due to 
the Euler-Maclaurin summation formula: 

N 



lim J- y l-^^^) 



k-h 
M 



k^hn 



\ M ) 



D 



N,M-^oo M 
N/M-^D 



For n > 2, we have 



k=l 
k^hM 



{k - huY 
which ends the proof. 



N 

E 

fe=i 

k^hu 



\f{x)-fixh)-{x-xh)f'{xh)\ 

{x - XhY 



(k - Hm^-^ (k - hM? 



< 



N 

E 

k=l 
k^hM 



l/( 



A' 

M' 



dx = 0. (A.3) 



/(^)l 



(k - Hm? 



(A.4) 

□ 
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Lemma A. 2. With the hypothesis of Lemma \A.l\ let no € N 6e such that sup^.g[o < no 

and \ f{xh)\ < nQ. Then 

\k-hM\>no 



Proof. Expanding the logarithms into their Taylor series we obtain 

f k \ _ fi 

{k - Km) 



N oo , -|xr+l fr(_k_\_ fr ( h-M \ 



k=l r=l 
\k-hM\>no 



The result of Lemma lA. II shows that only the terms corresponding to r = 1 give non- vanishing 
contributions. These are computed as a Riemann sum: 

■ M^r^ "O.JVW; AT A,,f A/f k~hM L T - ^ ' 



N,M ^oo ' N,M ^oo M ^ Jo X — Xh 

N/M^D N/M^D I, ^=1 

\k-hM\>no 



□ 



Corollary A.l. Let f and {Hm) satisfy the hypothesis of Lemma \A.l[ Let ^(A) be a strictly 
monotonous smooth function on M such that ^""^(0) = —q and ^(.D) = q. Define = 
(,~^{k/M), k G Z and p{X) = ^'(A). Then, independently whether the thermodynamic limit 
of Xhj^j belongs or not to the interval [—q,q], one has 



N 



1 



N,M^oolJ-k-hM + f{XhM) \J C(A)-C(Ah) 

N/M^D >^-^ \-q ) 

Proof. Assume that Hm C {1, • • • i and set / = / o Choose no > sup^.g[o,D] 1/(^)1 ^'^'^ 
decompose the original product into 

^ k-hM + f{\k) -TT k - hM + M) yr k - hM + M) 

^^k-hM + f(\h) k-hAr+f(!^) k-hM+Ji^V 

k=l ^ |fc-/iMl<no '''^'^ ^ J \ M I \k-hM\>no '''^^ ^ J K M ) 

The first product contains a fixed finite number of factors. Each of them goes to 1 in the 
limit considered. Hence, the first product goes to 1. The second can be calculated by using 
Lemma lA.2[ We obtain 

j.^ fr k-hM + fjXk) f fix) - fix,) J .^,0. 

where x, = lim hM/M. Equation (lA.Sp follows from (lA.lOp after the change of variable 

N,M^+oo 

A = ^^^(x). The case of Hm {1, ... , N} is proven along the same lines. □ 
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B The Fredholm determinant representations for the scalar 
products 

The proof of (14. 5|) , (14. 6|) is based on the following representation for a scalar product [Ml ESI HO] , 
valid whenever , . . . , /^^^^^ satisfy the system ()2.12p and z/i, . . . , z^at^ are generic complex 
numbers: 



n dijlla) 

(V'k(M)IV'(M)) = -j^ — detJ7.(W,MIM). (B.l) 

n sinhini^ - /i4) sinh(i/b - Ua) 

a>b 

The iV^ X Ni^ matrix defined as 



a=l 

Kd{iyj)t{uj,fiiJ sinh(/i<;^ - i/j + iC) , (B.2) 



a=l 

with 



-isinC , [ a (^) = sinh*^ (/i - </2) 

t(jU, z^j = . ^ . . . ^ . — ^ and < . (B.3) 

smh{^l - v) smh(^ -v-iQ i d (/i) = sinh^^ + iC/2) 



In order to obtain the scalar product (14.50 one should set here ^'j = \j for j = 1, . . . , A/^, where 
Aj are the Bethe roots describing the ground state. For the scalar product ()4.6p one should set 
in addition z^iv+i = —iC,/2. 

To obtain a Fredholm determinant representation for the scalar products we should present 
the original determinant of the matrix Vl^ in the following form 

AT 

dein,{{^^},{u]\{^l])\^^^^^=H{{^,},{\])dei(^^ JJ(A^._A,)-i). (B.4) 

a=l, 

Here r2(Aj, Afc|{/x}) is a new N x N matrix, H{{fj,}, {A}) an external coefficient. If we succeed 
to find a representation of the type ()B.4p . then we can replace the determinant of the N x N 
matrix by Fredholm determinant of an integral operator as 

det((^,fc + 0(A„Afc|{/i}) HiXj-Xay') = det (^I + —h{w,w'\{fi})ll{w - Xa^y (B.5) 

a = a=l 

Here the integral operator in the r.h.s. of (|B.5p acts on a counterclockwise oriented contour Tq 
surrounding the Fermi zone [—(7,(7]. This contour is such that it contains all the ground state 
roots Xj and no other singularity of the kernel Q{w,w'\{n}). 
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The proof of ()B.5p is quite obvious. Indeed, expanding the Fredholm determinant into the 
series of multiple integrals we see that each of these integrals reduces to the sum of the residues 
in the points Xj. Thus, the series of multiple integrals turns into the series of multiple sums. 
Then one can easily convince oneself that this series coincides with the expansion of the N x N 
determinant in the l.h.s. of (IB.5|) . 

Thus, our goal is to pass from the original representation ()B.ip to the form (IB.4p . The way 
we use here is based on the extraction of a Cauchy determinant from the original determinant 
defining the scalar product, an idea which was first introduced in [551 ES]- In the case of the 
scalar product (|4.5p this was done in [36j . The same method with minor modifications described 
below can be used for the proof of the second determinant identity (|4.6p . 

Following [SB], we consider the matrix ^}^{{z}, {i'}\{z}) for two sets of N,^ = + 1 generic 
parameters {z} and {z^}. In this case we write 



N+l 

A^ n sinhfzo 

det [ilf^l = —— — with Ajk 

AT+i detAT+i^ 



N+l 

Y\ sinh(zj 
a=l 



^ j coth(2;j — I'k) ioT k < N, 
^j \l forA; = iV+l. 

(B.6) 

The effect of multiplication by the matrix A is computed similarly to the method presented in 
[36] • We get 

det[n^{{z},{u}\{z})] = det ^ 



N+l 



N+l 



sinh (zk — Uj 



■ det Sik 
N+l ■' 



where 



c A 1/ n , na=i smh(i/fc-Za d 

Sjk = Sjky^iiyj\{z}) + . • w-y.ii^j\{y}) 

smh(i/fc - Ua) oyk 



k< N 



and we have set for arbitrary complex yi, . . . , yN+i 

N+l N+l 
yi^{v\{y}) = o-i^) n ^'^"^^(yk -V -K) + K-d{u) Jl sinh(yfc - v + iQ 

k=l k=l 



(B.7) 

(B.8) 
(B.9) 

(B.IO) 



Finally we reduce the size of det7v+i -S* by one by performing linear combinations of the lines 

SN+l,k 



det [Sjk] = Sn+i,n+i det 



N+l 



N 



Sjk ~ Sj^N^i 



SN-i 



1,N+1 



(B.ll) 



Thus, we arrive at the following representation for detAT-t-iO^ depending on generic complex 
zi,..., ZN+i and i^i, . . . , i^n+i- 



det [n^{{z}, {u}\{z})] = Sn+i,n+i det r 

N+l N+l [smh [Zk — Vj) 

where Sjk are given by (jB.Sp . (|B.9p . 



det 

N 



'>j,N+l 



Sn-\ 



l,k 



Sn+i,n+i 



(B.12) 
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Now we set Za = for a = 1, . . . , + 1, Vk = fo^^ A; = 1, . . . , and vn+i = — 

We assume that the parameters {^n^} satisfy the system ()2.1ip . while the parameters {\k} 
satisfy the system ()2.2p . Then we obtain for j < N, 

/"^.''^''^''^^ = Ksinh(A, - I) - sinh(A, + f ), 



116=1 sinh(Afc - \j - iC) 
dy,^ (lyj I {y}) /dyk 



n6ti sinh(Ab - Xj - <) 
and for j = N + 1, 



sinh(Aj — sinh(Aj + ^] 

(B.13) 



{y}={u} tanh(Afc - Xj + iC) tanh(Afc - Xj - iC) 



yn I {J^}) = a(-iC/2) sinh(-0 sinh(Ab - yj, 



N 



6=1 

N 



(B.14) 



dy. {uN+i I {y}) /dyk _ , = a(-</2) sinh(-iC) coth(Afc - |) TT '''^Hh " f 



It also follows from ()2.4p , (12.120 and the definition of the shift function ()2.13p that 

^ Vt sinh(/x^„ - w + iC) yr sinh(Aa - w - <) ^ ^2^*i?(,«) 
sinh(/x^^ — tt; — iC) sinh(Aa — w + i^) 

Substituting all these formulae into ()B.12p we arrive after simple algebra at the representation 
of the form ()B.4p . and thus, we obtain the Fredholm determinant representation (|4.6p . □ 
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